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ABSTRACT 

Context. High-mass stars, though few in number, play a major role in the interstellar energy budget and the shaping of the Galactic 
environment. However, despite this importance, the formation of high-mass stars is not well understood, due to their large distances, 
short time scales, and heavy extinction. 

Aims. Chemical composition of the massive cores forming high-mass stars can put some constrains on the time scale of the massive 
star formation: sulphur chemistry is of specific interest due to its rapid evolution in warm gas and because the abundance of sulphur 
bearing species increases significantly with the temperature. 

Methods. Two mid-infrared quiet and two brighter massive cores are observed in various transitions (E up up to 289 K) of CS, OCS, 
H 2 S, SO, S0 2 and of their 34 S isotopologues at mm wavelengths with the IRAM 30m and CSO telescopes. ID modeling of the dust 
continuum is used to derive the density and temperature laws, which are then applied in the RATRAN code to model the observed 
line emission, and to derive the relative abundances of the molecules. 

Results. All lines, except the highest energy SO2 transition, are detected. Infall (up to 2.9 kms~') may be detected towards 
the core W43MM1. The inferred mass rate is 5.8-9.4 10~ 2 M Q /yr. We propose an evolutionary sequence of our sources 
(W43MMl^IRAS18264-1152^IRAS05358+3543^IRAS18162-2048), based on the SED analysis. The analysis of the varia- 
tions in abundance ratios from source to source reveals that the SO and S0 2 relative abundances increase with time, while CS and 
OCS decrease. 

Conclusions. Molecular ratios, such as [OCS/H 2 S], [CS/H 2 S], [SO/OCS], [S0 2 /OCS], [CS/SO] and [S0 2 /SO] may be good indica- 
tors of evolution depending on layers probed by the observed molecular transitions. Observations of molecular emission from warmer 
layers, hence involving higher upper energy levels are mandatory to include. 
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1. Introduction 

OB stars are the main contributors to the evolution and energy 
budget of galaxies. Their formation, however, is not understood 
yet and the "class ical" scheme for low-mass star formation (see 
1 lAndre et"aT]|2000l and references therein) cannot be applied as 
such to OB stars. Indeed, young OB stars and protostars strongly 
interact with the surrounding massive clouds and cores leading 
to a complex and still not clearly defined sequence of objects 
from pre-stellar cores which are often believed to be hosted in 
the so-called IR dark clouds (IRDCs), to High-Mass Protostellar 
Obj ects (HMPOs), to hot cores and Ultra Co mpact HII regions 
(e.g. lBeufher et al.ll2007aUMenten et al]|2005l) . 

This tentative evolutionary scheme is however mostly based 
on observations of distant objects, for which, the physical scales 
are often at the order of 0.5 pc. In contrast, the well established 
protostellar evolution from Class to Class III Young Stellar 
Objects (YSOs) for low-mass stars revealed that the envelopes 
of individual protostars have typical sizes rangi ng from 0.02 to 
0.04 pc in clustered star-forming regions (see iMotte & Andrei 
I2OOH and references therein). HMPOs are therefore more proto- 
clusters than individual protostars. Recent interferometric obser- 
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vations have confirmed that they are fragmented and that they 
actually form cluster s of stars (e.g. IBeuther & Schilke l2004t 
IShepherd etall l2003h . More recently. lMotte et al.l d2007h have 
obtained a complete view of high-mass protostellar phases at the 
scale of a giant molecular cloud which is not too distant (Cygnus 
X at 1.7 kpc). This complex is rich enough to provide a first re- 
liable statistics (40 massive dense cores) at a physical scale of 
~ 0.1 pc. From this unbiased survey of the earliest phases of 
the high-mass star formation, it was deduced that massive pre- 
stellar dense cores are extremely rare, and that of the order of 
50 % of the massive dense cores are already forming high-mass 
stars while being still cold and not bright in the IR (IR-quiet 
dense cores). The remaining 50 % of dense massive cores can 
be referred as IR -bright dense cores. R eferring to the typical 
sizes discussed in IWilliams et al.l d2000l) . these 0.1 pc size ob- 
jects are dense cores, in contrast to clumps which are more 1 
pc size. On the other hand, these 0.1 pc dense cores are not 
necessary precursors of single stars and may form a group of 
stars, and would t herefore be con s idered as clumps in the pri- 
mary definition of IWilliams et ail d2000t>. We will here follow 
the scheme discussed in IMotte & Hennebelld (|2008) in which 
the study of massive dense cores at physical scales of ~ 0.1 
pc or less is required to really address the question of evolu- 
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tionary stages for high-mass star and cluster formation. The in- 
terest of studying the evolution of massive dense cores is also 
supported by the most recent theoretical scenarios, which pre- 
dict either a monolit hic collapse of turbulent gas at the scale of 
massive dense cores dTan & McKedl2002t iKrumholz & McKeel 
120081) . or a competitive accretion inside the gravitational po- 
tentia l of a cluster-forming massive dense core dBonnell & Batd 
120061) . Finally, it is worth noting that some direct observational 
evidence for the so-called global collapse has usually been ob- 
tained for massive dense cores at the typical ~ 0.1 pc scale (e.g. 
iMotte et al.ll2005l) . 

The chemical composition of the protostellar e nvelopes can 
be used to constrain the ages of low-mass protos t ars dMillar et al 



Table 1. Source Sample. The F21 value is the MSX value for 
the 21 jum flux, corrected for the distance of each source (F21 X 
id I \. If). 



19971 lHatchell et al.l 119981: iDotv etal.1 1200I IWakelam et al 



2004al) . The sulphur chemistry is of specific interest because 



of its rapid evolution in warm gas and since the abundances of 
sulphur bearing species increase significantly with temperature, 
both by ice evaporation and by shock interaction. Sulphur could 
thus act as a clock on time scales re l evant to the embedded phase 
of star formation (ICharnlevlll997l IWakelam et aill2004allbl) . A 
good sketch is now established for these objects allowing to date 
them within a "class" of objects. Here we propose some first con- 
straints on the time of evolution of the targeted massive cores; 
detailed chemical modelling will be presented in a forthcoming 
paper (Wakelam et al., in preparation). 

For high-mass protostars, the first studies using sulphur 
bearing mole c ules as e volutionary tracers are reported in 
lHatchell et all d!998l) and Ivan der Tak et al l d2003l) . We present 
here a complementary study of sources, whi ch are weaker at 
mid-in frared wavelengths. Also in contrast to Ivan der Tak et ail 
d2003l) who used c hemical m odels without the inclusion of 
CharnlevT 



atomic o xygen (e 



19971). we p r efer to use the recent 



Vastel et all (120021) . Furthermore, 
underlined that the sulphur may be 



results of lLis et al 2001t) and 
IWakelam et al.l J20041 " 
mainly in the atomic form in hot cores. 

To analyze these observat ions, we will apply a method sim- 
ilar to IWakelam et al.l (|2004a) who measured the age of the hot 
corinos of IRAS 16293-2422 with reasonable success. Our ob- 
servations of low- and high-energy transitions for each selected 
molecule enable us to probe the inner (> 100 K) and outer (< 
100 K) regions of the objects studied. Using transitions at sev- 
eral energie s, we first constrain, with a radiat ive transfer model 
(RATRAN. lHogerheiide & van der Takll2000l) . the sulphur bear- 
ing abundance profiles within the envelopes (this paper). In a 
forthcoming paper, the derived abundance profiles will be com- 
pared to a grid of chemical models (NAHOON) adapted to the 
physical structure of each source in order to constrain the age of 
each source. To cover the complete range of e nergy of the mas - 
sive envelopes (from 20 K to more than 100 K. lDotv et al.l2002l) . 
transitions with upper energies between 10 and 110 crrT 1 are re- 
quired. Indeed, even if observing emission lines covering a large 
energy range is mandatory to derive the excitation conditions of 
the molecules, the only way to probe the inner warmer regions 
is to observe high energy transitions which levels are not pop- 
ulated in the cold outer gas. Actually, probing the inner region 
of massive dense cores is the only way to chemically date their 
evolution from the molecular cloud. 



2. Description of the Source Sample 

The selected sources are two IR-quiet dense cores (MSX 
Flux2i um < 10 Jy at 1.7 kpc, following the defini- 
tion of IMotte et ail 120071) and two slightly brighter ones 
(IRAS05358+3543 and IRAS18162-2048), with bolometric lu- 





Source Coordinates" 




Properties' 1 


F 2 < 


Source 


RA (J2000.0) 


Dec 


Vlsr 




d 






( h m s ) 


C") 


(kms- 1 ) 


(is) 


(kpc) 


(Jy) 


IRAS05358+3543 


05 39 13.1 


+35 45 50 


-17.6 


6.3 X 10 3 


1.8 


18.9 


IRAS 18 162-2048 


18 19 12.1 


-20 47 31 


+ 12.1 


2.0 X 10 4 


1.9 


141.4 


IRAS 18264-1 152 


18 29 14.4 


-11 50 23 


+43.6 


1. X 10 4 


3.5 


4.9 


W43MM1 


18 47 47.0 


-01 54 28 


+98.8 


2.3 X 10 4 


5.5 


4.S 



" Coordi nates correspo n d to the peak of the mm c ontinu um emis- 
sion from iBeuther et al] d2002al) and IMotte et all d2003T) . respec- 
tively for IRAS05358+3543, IRAS18264-1152 and W43MM1. For 
IRA S18162-2048, coordinate s correspond to the CO emission peak 
from lBenedettini et alJd2004l) . 

b Luminosity, distance and velocities are fr omlSridharan et aljd2002h 
(IRAS05358+3543 and IRAS18264-1 152 ), IBenedettini et alTho04h 
(IRAS 18 162-2048) and lMotte et all f2003h (W43MM1). 



minosities 0.6-2.3 10 4 L at distances 1.8-5.5 kpc (see TableQ]) 
and sizes of ~0. 11-0. 13 pc (as seen in the beam of our obser- 
vations, see Sect. 16.11 ). These massive dense cores are expected 
to be fragmented at small scales and are meant to represent the 
earliest phases of the formation of clusters. The expected mul- 
tiplicity can however be neglected to first order, to follow the 
average chemical evolution of their surrounding dense cores. In 
fact, a small fraction of the total mass of the cores is in the in- 
dividual fragments, and usually the large scale average density 
profile of the cores applies well down to the smallest observed 
structures. This suggests that the cores could be well represented 
to first order as spherically symmetric structures from a few 0.1 
pc down to ~ 0.01 pc. IRAS05358+3543 has for instance bee n 
studied with millimeter interferometers (see lLeurini et al.l2007l) . 
It appears that the ~ 0.3 pc size core of 613 M splits into two 
main objects (mmla and mmlb) of 1 and 0.6 M inside regions 
of 0.006 pc diameter, separated only by ~ 0.008 pc. A mass 
of 1.6M in a region ~ 0.01 pc comprising mmla and mmlb, 
inside a 0.3 pc size core of 613 M corresponds to an average ra- 
dial profile of the density of 25 which is actually clo s e to the 
r~ L4 profile derived from large scale by IBeuther et al.l d2002bl) 
and which we use to describe the global ID distribution of mat- 
ter in the core. Two other millimeter condensations mm2 and 
mm3 are also detected with the interferometer. They could even- 
tually confuse the results for molecular tracers which can trace 
such high density regions. 

A large number of continuum and molecular line (CO, SiO, 
CH3OH, H2 CO etc) data towards th ese sources have been pub- 
lished {e.g. IBenedettini et al.l 120041) . The physical conditions 
have thus already bee n determined in some of the envelopes {e.g. 
IBeuther et alj l2002bT) and we use the published data together 
with our observations. We limit our study to massive dense cores 
to simplify this first step of the analysis for two reasons. First, 
different physical conditions have been observed in these ob- 
jects, but they still have the same structure in the sense that there 
is no strong UV field or ultra-compact HII regions, which will 
change drastically the chemistry with time. The second reason 
is that these cores could be for massive dense cores the equiva- 
lent of the Class YSOs in the general scenario of gas and du st 
evolution, for which sulphur chemistry (Wa kelam et al.l2004allbh 
has already been studied. Finally, the sources studied here are 
part of the sample that will be studied within the key program 
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Fig. 1. Spectral energy distributions obtained from ID model 
overlaid on observed fluxes of each source. Peak fluxes in the 
millimeter range have been adjusted to fit radial extension of the 
model. In each caption are given the temperatures of the black 
body components (top right), the maximum flux F max and the 
corresponding wavelength of the SED, the flux F\i at 12//m, and 
the contribution of the hot part (F35, integrated flux for A <35 
yum) to the total integrated flux F tota i (top left). 



W/5/£] with the HIFI instrument aboard the ESA Herschel Space 
Observatory. 

Both IRAS05358+354 3 and IRAS18264-1 152 are part of 
the lSridharan et alj d2002l) sample. IRAS05358+3543 is a rela- 



tively low luminosity (L ~ 6.3 10 L Q ) and nearby (1.8 kpc) 
massive dense core, composed of three main sources (mml, 
mm2 and mm3; tota l mass around 613 M ), all within 4"x9" 
(iLeurini et al.l 120071) . hence within the telescope beam. The 
source exhibits a hyperc ompact HII region wh ich may evolve 
into an UC-HII region. IBeuther et ail d2007bl) show that the 
source actually splits into four individual continuum compo- 
nents, two of which are part of a proto- binary system, mml (dy - 
namical age ~ 3.6 10 4 yr). According to lLongmore et al .| ( 20061). 
the mai n mid-IR source coincides with mm la of IBeuther et aD 
(I2007bl) . but the second mid-IR source does not coincide with the 
mm source mmlb, hence suggesting even more sub-structure. 
The source exhibits a class II methanol maser indicating that it is 
a very young object. The source mm3 is the coldest and youngest 
(Tdust ^ 20 K). Several molecular outfl ows are observed (tw o 
from mml, outflow rate of 6 x 10 4 M G ; IBeuther et ai]l2002al) . 
St rong line em i ssion asymmetries in HCO + have been observed 
bv lFulleretai]d2005l) . 

IRAS 18264-1 152 is a massive dense core at a distance of 
3.5 kpc. Its luminosity is estimated to be up to 10 4 L . A 
massive molecular o utflow is centered on the mm dust peak 
dBeufher et alJl2002al) . Continuum observations with Plateau de 
Bure dOiu et al.ll2007l) reveal that the source is double; the mass 
of the dust-gas core is estimated to be 570 M , while the outflow 
ma ss is 20.2 M n (dynamical age of 0.5 10 4 yr). 

iBenedettini et al.l d2004l) performed multi-line observations 
of IRAS 18 162-2048, a luminous infrared source (L= 2 10 4 L Q ; 
lYamashita et alj Il989h at a kinematic distance of 1.9 kpc 
dKurtzet al.lll994l)~ This object is associated with the IRAS 
source and is actually an IR-bright dense core. A very power- 
ful outflow (M ~570 M ), one of the most massive known, was 
observed (dynamical age ~ 10 6 y r). The source has bee n resolved 
into a cluster of several objects dStecklum et al.|[l997l) . 

W43 is a massive star-forming region, located at 5.5 kpc, 
which harbors a g iant HII region and was studied in detail by 
iMotte et ail d2003l) . W43MM1 is the most massive dense core 
of the region and does not coincide with any known IR source. 
The dust temperature is estimated to be 19 K and its luminosity 
2.3 10 4 L Q . W43MM1 harbors a methanol and a water maser, 
but does not show any infrared or centimeter emission. An infall 
has been detected by Motte et al. (private communication). 

Continuum emission from our sources is quite well 
determined with well sampled observ ations from IRAS , 
MSX, SCUBA/JCMT or SHARC/CSO dHunter et al.l 120001: 
Uenness et all ["1995; Mi nier et all 120051) . SMA dBeufher etal 



2007bHSu etai]|2004t) and MAMBO/IRAM30m dBeufher et al 

2048 is particularly well con- 
complete ISO-SWS observation. Fewer 



http://www.strw.leidenuniv.nl/WISH/, PI: E. van Dishoeck 



2002bD . The SED forIRAS18162- 
strained, thanks to a 
continuum observations are available for the source W43MM1. 

From the fitted SEDs (see Sect. 16.11) . shown in Fig. Q] we 
make a rough evolutionary classification of our 4 objects, using 
the following parameters: 

- wavelength and flux of the maximum continuum emission, 

- flux at 12 fim, 

- contribution of the hot part (A <35 //m) to the total integrated 
flux, 

- temperature of the black-body components. 

We assume that the less evolved the source is, the coldest it is, 
hence the SED peaks at longer wavelength with weaker flux. As 
the massive core evolves, it becomes warmer hence heating the 
dust wherein it is embedded: as a consequence the contribution 
of the flux at shorter wavelength (F35, the integrated flux for 
A <35 jum) increases. 



4 



F. Herpin et al.: S-bearing molecules in Massive Dense Cores 



Applying these criteria, W43MM1 appears to be the 
youngest object. A sequence of evolution W43MM1 — > 
IRAS18264-1152 -> IRAS05358+3543 -> IRAS 18 162-2048 
is proposed and adopted for the following discussion. This evo- 
lutionary sequence is also suggested by the dynamical ages de- 
rived by the authors cited above. 



3. Observations 

The observations presented here were performed during two 
sessions. The first one in September 2005, using the IRAM- 
30m antenna^. The second one in May 2006, using the Caltech 
Submillimeter ObservatorjQ, on the summit of Mauna Kea 
(Hawaii). 

With the 30m telescope, all sources were observed in the ro- 
tational transitions lines of CS 7=5-4, SO 7=34 - 23, 7=56 ~ 
4s, 7=65 - 5 4 , S0 2 7=5 1<5 - 4o, 4 , 7=10 ,io - 9 h9 , 7=11 U1 - 
10 ,io, 7=14 ,i4 - 13 U 3, OCS 7=8 - 7, 7=13 - 12, 7=19 - 18, 
H2S7=lio - Ion J-^-iq - 2\ 1, and in the isotopic lines of 
34 SO 7=3 4 -23, S0 2 7=5i 5 -4 4 , OC 34 S 7=8 - 7, H 34 S J=h - 
l ,i (see Table Only IRAS05358+3543, IRAS 1 8264- 1 i 52 
and W43MM1 were observed in the CS and C 34 S 7=3 - 2 
lines. Observations were made using multiple receivers (either 
simultaneously A100, B100, A230, B230 or CI 50, D150, C270, 
D270) coupled to the high-resolution VESPA backend (resolu- 
tion of 80 kHz) in frequency switching model (frequency throw 
of 7.9 MHz). The sky conditions were reasonable (r? m ~0. 1-0.5 
at 225 GHz, see Table [2] for the T sys ). The telescope pointing 
(better than 2") and focus have been set on suitable planets (i.e. 
Jupiter) or standard calibration sources. 

Complementary observations have been performed with 
the CSO telescope for lines with frequencies higher than 
300 GHz (in order to better cover the energy ladder). As 
IRAS05358+3543 was not visible at the CSO during our run, 
only the other three sources were observed in lines of CS 7=7-6, 
S0 7=8s - 77, SO2 7=18oi8 - 17i,n, and in corresponding iso- 
topic lines of C 34 S, 34 SO, 34 S0 2 . The receiver is a helium-cooled 
SIS mixer operating in double-sideband mode designed to fully 
cover the 280-420 GHz window. The backend consists of a 1024 
channel acousto-optical spectrometer covering a bandwidth of 
50 MHz, providing a velocity resolution of Av 0.086 km s at 
345 GHz. Observations were performed under good but unsta- 
ble weather conditions (if" ~ 0.06-0. 15 at 225 GHz, see Tabled 
for the T sys ). The pointing accuracy was of the order of 4-5". The 
CS 7=7 - 6 line was blended with an H 2 CO line from the im- 
age sideband, therefore we applied a frequency shift to the local 
oscillator frequency to separate the two lines, but the procedure 
did not work well for IRAS 18264-1 152 and W43MM1 (hence 
the CS spectra are contaminated). 

Data reduction was performed using the CLAS S software 
from the GILD AS suite dGuilloteau & Lucas! I20001) . We found 
and eliminated unusable data, subtracted a baseline, then spec- 
tra at the same position were summed and finally antenna tem- 
perature T* was converted into T m \, (using the i] m b values from 
Tab. |2|. Special attention was paid to the frequency switching 
observing method. 



2 http://iram.fr/IRAMES/index.htm 

3 http://www.cso.caltech.edu. The Caltech Submillimeter 
Observatory is supported by the National Science Foundation un- 
der award AST-0540882. 
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V LSR( km / s ) 

Fig. 2. Normalized emission of CS, H 2 S, OCS, SO, S0 2 (and 
isotopic species) lines from W43MM1. Spectra velocity reso- 
lutions are 0.10-0.19 kms -1 . The dashed line shows the source 
LSR velocity. 



4. Results 

We present here the result of the observations for each source. 
Some spectra of the detected lines are displayed in Fig.[2f|5](the 
remaining spectra are available as online material). Some tran- 
sitions ( 34 S0 8 8 - 7 7 , S0 2 28 4 ,24 - 283,25 and 34 S0 2 5 lj5 - 4 , 4 ) 
are not detected in any of our sources. The line parameters are 
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Table 2. List of observational parameters for the IRAM 30m and CSO telescopes (observed species, energy level transitions, 
upper energy, line emission frequencies, half power beam width, instrument, main beam efficiency 7] m b, receiver name, velocity 
resolution 6u and system temperature r sys . 



Species 


Transition 


up 


Frequency 


HPBW 


Instrument" 




Receiver 


8v 


f 






[cm- 1 ] 


[GHz] 


["] 








[m-s- 1 ] 


[K] 


CS 


3-2 


10 


146.969 


17 


30m 


0.69 


D150 


160 


315-350 




5-4 


24 


244.936 


10 


30m 


0.48 


A230 


96 


530-860 




7-6 


46 


342.883 


24 


CSO 


0.75 


345 


43 


1210-1340 


C J4 S 


3-2 


10 


144.617 


17 


30m 


0.69 


D150 


162 


270-390 




7-6 


46 


337.396 


25 


CSO 


0.75 


345 


43 


1260-1760 


so 


3 4 -2 3 


11 


138.178 


18 


30m 


0.70 


D150 


170 


300-330 




5 6 -4 5 


24 


219.949 


11 


30m 


0.55 


B230 


110 


590-930 




65 — 54 


35 


251.825 


10 


30m 


0.48 


D270 


190 


750-830 




8 8 - 7 7 


61 


344.308 


24 


CSO 


0.75 


345 


43 


1020-1050 


,4 SO 


3 4 -2 3 


11 


135.775 


18 


30m 


0.70 


D150 


170 


290-350 




8 8 -7 7 


61 


337.582 


25 


CSO 


0.75 


345 


43 


1020-1060 


so 2 


5 1 5 - 4o4 


11 


135.696 


18 


30m 


0.70 


D150 


170 


270-400 




lOo'io - 9 1>9 


34 


160.827 


15 


30m 


0.66 


C150 


150 


350-520 




1 1 1.11 - lOo.io 


42 


221.965 


11 


30m 


0.55 


B230 


100 


530-550 




14o,14 _ 13l,13 


65 


244.254 


10 


30m 


0.48 


A230 


190 


830-1430 




18o,18 - 17l,17 


105 


321.330 


26 


CSO 


0.75 


345 


45 


770-820 




284 -2 4 - 283,95 


289 


267.719 


9 


30m 


0.46 


D270 


170 


470-510 


,4 so 2 


5l,5-4o,4 


11 


133.471 


18 


30m 


0.70 


D150 


170 


200-220 


ocs 


8-7 


14 


97.301 


25 


30m 


0.77 


A100 


240 


140-160 




13-12 


37 


158.107 


15 


30m 


0.66 


C150 


150 


340-360 




19-18 


77 


231.060 


11 


30m 


0.52 


A230 


100 


450-580 


OC J4 S 


8-7 


14 


94.922 


25 


30m 


0.77 


B100 


240 


140-160 


H 2 S 


1 1,0 - lo.i 


19 


168.762 


14 


30m 


0.70 


D150 


140 


610-800 




2 2 ,o _ 2^1 


58 


216.710 


11 


30m 


0.57 


B230 


110 


410-600 


H 2 34 S 


li.o _ lo.i 


19 


167.910 


14 


30m 


0.70 


D150 


140 


630-670 



" Conversional factor is S /T„ lb = 4.95 Jy/K for IRAM 30m telescope and S /T mb = 74 Jy/K for CSO telescope. 



given in Table [6H9] (only available in electronic form). For all 
observations, contamination from lines in the image sideband of 
the receiver has been checked. 

4.7. W43MM1 

The CS line emission from this source is characterized by what 
looks like a strong self-absorption (at ~ -0.2 kins -1 from the 
source velocity) with a stronger blue component compared to 
the red one. Indeed the self-absorption in the 3-2 line goes 
below zero kelvins, splitting the line profile into two well sep- 
arated components. The 2mm continuum flux is actually high 
(0.8 Jy), which converts into a brightness temperature of 5.5 K, 
and can therefore explain the observed dip. Similar, but weaker, 
self-absorption is observed in the 5-4 and 7-6 transitions. 
However, the isotopic lines are only weakly asymmetric, red- 
shifted with respect to the systemic velocity. The derived opaci- 
ties (from our model, see Sect. 16.2b are respectively 0.4 and 0.6 
for the C 34 S 3-2 and 7-6 line emissions, indicating optically 
thick lines. All lines are very broad, up to 15 km s" 1 . Note that 
HCO 4 and H 13 CO + 3-2 emission lines observed bv lMotte et"ai] 
(120031) have exactly the same profiles as the CS and C 34 S emis- 
sion lines, respectively. The self-absorption is present in the FLS 
line emissions too (the isotopic line is not detected), even pro- 
ducing for the 1 1,0 — lo,i line a negative dip, below the continuum 
level, at the same velocity as CS. The 2n^-2\^\ line emission has 
a different profile with more redshifted emission. 

The SO 34-23, 56-45 and the isotopic 34-23 line profiles are 
also weakly asymmetric, with a more pronounced blueshifted 



emission. The 65 - 54 and 8s - I7 lines exhibit, on the contrary, 
a stronger redshifted component. The SChSi.s - 4o,4 line is also 
blue-asymmetric, but all other SO2 lines are more or less sym- 
metrical and gaussian. 

All OCS lines exhibit a gaussian profile, slightly red-shifted 
by ~ 0.5 kms -1 from the source velocity, with a red wing. The 
isotopic OC 34 S line profile is complicated with several compo- 
nents in emission at 91.8, 94.0 and 95.5 kms -1 (see Fig. [2j. No 
identification with known species was possible (we also checked 
the frequencies in the image band) and instrumental effects have 
been ruled out. This weak component is not seen in the corre- 
sponding OCS line profile. Since W43MM1 is a very rich molec- 
ular source (see Sect. 17. II ). this could well be real lines. 

4.2. IRAS1 8264-11 52 

All line emission profiles are well peaked on the source veloc- 
ity and have more or less the same shape, with a blue wing. 
Nonetheless, a distinct weak emission at 48-50 kms -1 is ob- 
served for all transitions of all species. These emissions cannot 
be due to a contamination from any OFF position as observations 
have been performed in frequency-switch mode. None of these 
features correspond to a known molecular line. We propose that 
this secondary emission comes from a different location than the 
main one, either from another source within the beam, or from 
a knot in the outflow. Actually, an SiO out flow component a t 48 
kms"' 1 , offset by 20" has been observed by lOiu et al.1 (120071) . 

All CS lines are detected, except the 7-6 transition that is 
blended with an H2CO line from the image band, hence com- 
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Fig. 3. Normalized emission of CS, H 2 S, OCS, SO, S0 2 (and 
isotopic species) lines from IRAS18264-1 152. Spectra veloc- 
ity resolutions are 0.10-0.19 kms" 1 . The dashed line shows the 
source LSR velocity. 



pletely contaminating the profile. The C S 7-6 line is tenta- 
tively detected. The H 2 S lines exhibit a less peaked profile with 
a blue wing too. All the SO lines are detected too, but the 65 - 54 
line emission exhibits a strong and very broad red wing which 
looks suspicious. We do not consider this line for further analy- 
sis. The SO2 line profiles are more flattened. Actually, an uniden- 
tified line at 47.8 kms 1 is observed in the SO2 5i,5-4o,4 spectra, 




-25 -20 -15 -10 

V LSR( km / s ) 



Fig. 4. Normalized emission of CS, H 2 S, OCS, SO, S0 2 (and 
isotopic species) lines from IRAS05358+3543. Spectra veloc- 
ity resolutions are 0.10-0.19 kms 1 . The dashed line shows the 
source LSR velocity. 



and could come from the reference OFF spectra, just as the sus- 
picious 65 - 54 wing. The OCS lines exhibit similar profiles, with 
a blue asymmetry. 
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Fig. 5. Normalized emission of CS, H 2 S, OCS, SO, S0 2 (and 
isotopic species) lines from IRAS18162-2048. Spectra veloc- 
ity resolutions are 0.10-0.19 kms" 1 . The dashed line shows the 
source LSR velocity. 



4.3. IRAS05358+3543 

The CS line emission is strong in both observed transitions (the 
3-2 line emission being the strongest and broadest one) with 
a strong dip at the source velocity, even for the isotopic line. 
The isotopic line emission here is optically thin (around 0. 1 ac- 
cording to our model), therefore this dip is likely owing to two 



velocity components rather than to self-absorption. All profiles 
are asymmetrical with a stronger red component, the blue one 
being almost completely extinguished (nearly in absorption) for 
the CS 5 - 4 line emission. The peak-to-peak velocity separation 
is 3.0, 3.5 and4.3 kmr 1 , respectively for the C 34 S 3-2, CS 3-2 
and 5-4 lines. 

The HtS li o - lo i line emission also exhibits a strong dip 
at the source velocity. The separation between the two emission 
peaks is 3.6 kms -1 . The same shape is also seen for the 22$—2\ t \ 
line, where no blue peak is seen, the red component being split in 
two peaks. The isotopic line might be tentatively detected, with 
a very strong asymmetry. 

Profiles of the SO lines generally exhibit a dip at the source 
velocity. The peak-to-peak separation is 3 kms"' for the 34 - 
23 line. The isotopic 34 - 23 line shows a red-shifted (by more 
than 1 km s" 1 ) double-horn profile. 34 SO 34 - 23 line emission is 
marginally detected, suffering from obvious baseline problems. 
We checked, without success, if the double-peak profile could be 
due to other species. 

The SO2 emission is detected in all transitions (except the 
11 i,u - 10o,io), but the S /N is quite low (5-9). All the lines 
are asymmetric, nearly without any red component. The 140,14 - 
131,13 emission is the strongest. 

A dip is seen in all OCS lines too, most prominently in the 
13-12 transition. The 19-18 line is the strongest. The OCS 
13-12 line is only tentatively detected with a complex profile. 

A second emission component is seen in the C 34 S 3-2, OCS 
8-7, SO 3 4 - 2 3 , S0 2 5i, 5 - 4 , 4 spectra at -13 kms" 1 . As this 
emission is also present in the wing of the CS 3-2 and H2S 1- 
1 line emission (observed with a smaller telescope beam), we 
believe that this emission originates from the source mm lb as 
observed in the C 34 S 5-4 line emission bv lLeurini et al.l d2007l) 
at the same velocity. 



4.4. IRAS1 81 62-2048 

The CS line shapes are similar for the different transitions ob- 
served, with the 5-4 line being substantially stronger and 
broader, exhibiting a red wing. 

The H2S emission is slightly broader than CS, with more 
complex profiles: the 1 1,0 - lo.i emission profile is made of two 
peaks, possibly revealing self-absorption at a redshifted velocity 
(by 0.5 km s" 1 with respect to the source systemic velocity). The 
22,o - 2 1,1 line profile is slightly asymmetrical, with a red wing. 
The optically thin isotopic line emission exhibits a symmetrical 
profile centered on the source velocity, with a narrower profile 
than for the H2S emission. 

The SO lines are strong and easily detected, but the isotopic 
line is not detected. The SO2 lines are well detected, even the 
18o, is - 17i t n transition, the lli.n - 1 Oo, 1 being the strongest 
one. The line profiles show a red wing. 

All the OCS lines are completely asymmetrical with a very 
weak blue component, indicating a self-absorption at the source 
velocity. 

For that source too, emission a few kms" 1 away from the 
main emission is detected (at 16-17 kms" 1 ) for the lines ob- 
served with a beam larger than 18", likely revealing the presence 
of a second source within the beam. 
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Table 3. Asymmetry parameter dMardones et al.|[l997l) calcu- 
lated for the detected line emissions. Negative values of 6v cor- 
respond to blueshifted emission while positive values are red- 
shifted. 



Species 


Line^ 




Source 










W43MM1 


18264 


05358 


18162 


CS 


3-2 


-0.69 


0.05 


0.37 






5-4 


-0.41 


0.16 


0.42 


0.03 




7-6 


-0.40 







-0.01 


C J4 S 


3-2 


0.05 


-0.02 


0.27 






7-6 


0.18 


0.20 




0.13 


so 


3 4 -2 3 


-0.10 


0.00 


0.22 


0.03 




5 6 -4 5 


-0.06 


0.05 


0.29 


0.06 




5 - 3 4 


n 1 1 


c 


U.zo 


-U.Uo 




8g - 7 7 


0.11 


o!o4 




0.12 


so 2 


5 1,5 - 4(1,4 


-0.11 


0.08 


0.17 


0.03 




IOq.io - 9,,, 


0.05 


-0.02 


0.27 


0.05 




11 l.ii _ 10o,io 




-0.03 




0.44 




14o,14 - 13l,l3 




-0.12 


0.38 


0.00 




I80.I8 - 17l,17 


0.10 






0.19 


ocs 


8-7 


0.11 


0.06 


0.26 


0.53 




13 - 12 


0.03 


-0.06 


0.24 


0.52 




19-18 


0.10 


0.09 


0.17 


1.08 


H 2 S 


1 1,0 - lo.i 


-0.76 


0.02 


0.39 


-0.37 




22,0 _ 2i,i 


0.06 


0.08 


0.42 


0.16 



a Reference lines are optically thin 34 SO and H 2 34 S line emissions 
* doubtful line 
doubtful line 



5. Analysis 

5.1. Line Asymmetries 

Analysis of the line profiles can reveal the presence of outflows, 
infall or even rot ation because of their known signatures (see 
iFuller et al.l I2005 , and references therein): outflow or rotation 
give rise to both red and blue asymmetric lines, while infall pro- 
duces only blue ones. Indeed, the profiles of optically thick lines 
from infalling material have stronger blueshifted emission than 
redshifted one. However, it is important to stress that outflow or 
rotation could also produce a blue asymmetric line profile along 
a particular line of sight to a source. 

To quantify the asymmetry of a line, we use lMardones et ail 
( 1 1997b criteria: we calculate the asymmetry parameter 



6v = 



^thick — fthin 



Al- 



l/2 



where i> t hi c k and i> t hi n are the line peak velocity respectively of an 
optically thick and thin tracer, and Avi/2 the line width (FWHM) 
of the optically thin line. 

Because not all the same lines are detected or observed in all 
sources, the difficulty is to find a common strong optically thin 
line, generally from isotopic species. We use 34 S034 - 23 for 
IRAS05358+3543 (v thin = -17.0 kms -1 , Au 1/2 = 4.4 10ns- 1 ), 
IRAS 18264- 1152 (u^ = 43.6 kms -1 , Au 1/2 = 4.0 kms -1 ), 
and W43MM1 (u lhin = 99.2 kms -1 , Av l/2 = 5.5 kms -1 ), and 
H 2 34 Sli, - l .i for IRAS18162-2048 (v thin = 12.1 kms -1 , 
Ai>i/2 = 1.5 kms -1 ). Results are given in Table[3] 

Negative values of 5v correspond to blueshifted emis- 
sion while positive valu es to redshifted one. Following 
iMardones et alj d 19971) and IFuller et alj d2005l) . we adopt a cri- 
teria of \6v\ > 0.25 to indicate that a line profile is asymmet- 
ric. Main conclusions are that the source IRAS 18264-1 152 has 



quite symmetric lines, except for the C 34 S 7-6 line (close to 
the 0.25 value), while for IRAS18162-2048 asymmetric lines 
are observed in OCS, in the H 2 S li o - lo, 1 (redshifted) and in 
the S0 2 11-10 transitions (blueshifted). Very strong asymme- 
tries are observed for the two other sources: most of the lines 
from IRAS05358+3543 are strongly redshifted (6v up to 0.42); 
W43MM1 shows clearly blueshifted lines in CS (up to -0.69) 
and H 2 S (up to -0.76), probably indicating the presence of in- 
fall towards this source. Nevertheless, we stress again that such 
an emission profile can be sometimes produced by an outflow 
too, e.g. CO 1 - line profile towards IRAS 18 162-2048 by 
iBenedettini et a l. (2004). Hence, this option cannot be elimi- 
nated without further investigations. 

5.2. Possible Infall in W43MM1 

Assuming the observed CS profile towards W43MM1 is due 
to infall, we estimat e the infall v e locity using the approximate 
method described in iMvers et~ai] dl996l) for the optically thick 
lines with a strong blue asymmetry: 



Vied - fblu 



■In 



1 + eT BD /T D 
1 + eT RD /T D 



(2) 



where To is the brightness temperature of the dip (assumed op- 
tically thick), Tbd and Trd the height of respectively the blue 
and red peak above the dip, and !%i ue and u re d respectively the 
velocity of the blue and red peak. The velocity dispersion <x is 
obtained from the FWHM of an optically thin line. 

Adopting the 34 SO 34 - 23 line as an optically thin tracer, the 
velocity dispersion in the circumstellar material is estimated to 
be 5.5 km s -1 and the source velocity to be 99.2 km s -1 . The de- 
rived infall velocities v,„ are 2.1, 1.8 and 2.9 kms~', calculated 
respectively for the CS 3 - 2, 5 - 4 and H 2 S 1 1,0 - lo, 1 line emis- 
sions. Using the following formula (where m is the mean mass 
of the molecule and r oM the inner radius of the most outer layer): 



dM 
dt 



= Anr olll 2 mnVi, 



the inferred kinematic mass 
3.4 10 17 cm and a density n of 2.5 
5.8-9.4 10~ 2 M /yr. 



infall rate 
10 6 



(for 

-3 



cm 



(3) 

a size r out of 
cf. Table 4) is 



5.3. Outflows in the other sources or multiplicity effect ? 

The SO 3 4 -2 3 lines towards IRAS18162-2048 can be fitted with 
(1) two Gaussians of FWHM line widths of 1.2 and 2.6 kms" 1 (see 
Fig. [6]). The CS 7-6 line emission can be fitted by two Gaussians 
too: one of 1 kms 1 broad, the other one of 2.8 kms 1 . The 
optically thin H 2 34 S line can be fitted by two Gaussians of 
FHWM 1.2 and 3.2 kms 1 , comparable to the components de- 
rived from the other lines. Moreover, it is interesting to note 
that the two gaussian components are centered at velocities 
(11.9 and 12.3 kms -1 ) comparable to the velocities (11.9 and 
12.8 kms 1 ) of the two major outflow components seen in CO , 
13 CO and C ls O 1 - line emissions bv lBenedettini et all d2004l) . 
For IRAS 1 8 1 264, the CS 3 - 2 and SO 3 4 - 2 3 line emissions can 
be fitted with two Gaussians of two different FWHM: 1.7 and 
3.7 kms -1 . These two-components reveal the presence of gas 
flows, possibly an outflow towards objects IRAS18162-2048 
and IRAS18264-1152. 

Line emission profiles seen in the source IRAS05358+3543 
exhibit a strong dip (especially for CS), even for the optically 
thin isotopic species, which is hence not a self-absorption but 
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Fig. 6. Gaussian fit of the SO 34 - 2j line emission from 
IRAS05358+3543, IRAS18162-2048 and IRAS 18264- 1152. 
The observed spectra are in black, total fit in red, individual com- 
ponents in blue and green. 



rather a the signature of two well separated velocity compo- 
nents with FWHM around 2.6 and 2.9 kms~' (derived from 
CS 3 - 2 and SO 34 - 23 lines emission fitting, see Fig. [ 6|. This 
cou ld be due to the outfl ow detected by IBeuther et al.l l2002al) 
and IBeuther et al.l d2002cl) or to the multiplicity of the source it- 
self. IRAS05358+3 543 is indeed composed of three mm sources 
dLeurini et alj|2007l) . all within our beam, mml being the most 
massive, the other two being less evolved (mm3 is probably a 
starless massive core). Actually, the comp onent at -16 kins' 1 
is at the same velocity as the source mm2 dLeurini et al. I I2007I) . 
Hence the profile may result from the source multiplicity. 




Fig. 7. Plot of the density (top) and temperature (bottom) versus 
the distance to the center for the 4 objects. 



parameters, we restrict our modeling to a simple ID, radial de- 
scription of the sources. The assumed spherical symmetry of the 
sources is acceptable if we restrict our study to the molecular 
line emission modelin g coming from cold gas (T < 100 K) as 
[Marseille et al.ld2008l) showed it for the massive dense core case. 



6. Modelling the continuum and the molecular 
emission of the sources 



We use the modeling metho d described inlMarseille et al.l ()2008l) 
and originally developed bv lHogerheiide & van der Tak d2000l) . 
It consists, in a first step, of constraining the physical structure 
of the source using a modeling of the SED, and in a second step 
of applying the obtained physical model to derive abundances 
of the observed species. In order to reduce the number of free 



6.1. SED distribution 



We use the radiative transfer code MC3D dWolf et al.ll999l) in its 
ID version to model the SED. The model parameters are taken 
from literature (distance, d, power law index, p, for density dis- 
tribution, heating source temperature, T», and source sizes, r out , 
the total luminosity, L). The parameter, p, which defines the den- 
sity distrib ution through a power law i n the form p ( r) oc r p , was 
derived by IBeuther et all d2002d) and iMotte et alj d2003l) from 
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v lsr (km.s ) v lsr (km.s ) 

Fig. 8. Fit of the CS lines emission from IRAS 1 8264- 1 152. The 
observed spectra are in black, fit in red. 



fits to the millimeter continuum emission maps. Concerning the 
source size, we adopt systematically the size of the continuum 
mm- wave emission, deconvolved by the beam. These values are 
summarized in Table [4] The actual value used for the inner ra- 
dius of the physical structure has virtually no observable effect 
on the results, we therefore adopt for the inner radius (r su b) the 
radius of dust sublimation taken at 1500 K. 

The fit of the millimeter and submillimeter SED of each 
source allows deriving the best temperature and distributions by 
iterating on the total mass which is controlled by the value of the 
density «o at reference radius ro taken equal to 100 AU, and on 
the average temperature. The values obtained, with mean densi- 
ties and temperatures of the sources, are reported in Table [4] and 
shown in Fig. [7j 

6.2. Modeling of the molecular emission 

The physical description of the sources obtained from the 
SED modelling is then used to model the molecular emis- 
si on. We use the RATRAN radiativ e transfer code developed 
bv lHogerheiide & van derTakl(l2000l) . We build data cubes with 
a velocity resolution equivalent to the observations, and with a 
high spatial resolution (0.5") which are convolved with the beam 
sizes of the observations for a direct comparison with the ob- 
served spectra. We adopt constant abundances for the RATRAN 
modeling. The abundances of all species are derived indepen- 
dently for each observed transitions of each species. Then the 
different abundances obtained for each particular species are in- 
terpreted as due to change in abundance as a function of the ex- 
citation temperature and therefore as a function of depth into the 
sources. 



Table 4. Summary of results from dust continuum emission 
modeling (above mid-line) and from molecular line emission 
modeling (below mid-line). Different turbulent (vj) or infall 
(tfinfaii) velocities are used to reproduce the different line emis- 
sions from each source. 



Source 


W43MM1" 


18264 


05358 


18162 


d (kpc) 


5.5 


3.5 


1.8 


1.9 


L(1O 4 L ) 


2.3 


1.4 


0.7 


2.9 


1 * (1U Iv) 


3.0 


i.y 


2.0 


3.0 


r out (10 4 AU) 


2.5 


2.7 


2.1 


2.6 


pi 


-2.0 


-1.5 


-1.4 


-2.1 


<(10 4 ) 


40 


1.0 


0.3 


6.1 


''sub ( AU ) 


4.9 


26.5 


20.1 


27.2 




-0.54 


-0.61 


-0.62 


-0.60 




3.55 


4.05 


3.99 


4.04 


<r> (K) 


20.0 


32.6 


31.1 


35.6 


Tout (K) 


14.2 


24.6 


24.1 


27.2 


(n) (10 6 cm" 3 ) 


11.8 


2.6 


0.2 


1 


M(M Q ) 


4100 


1200 


400 


570 


v, (km/s) 


0.8-1.5 


0.7-1.3 


0.7-2.6 


0.5-1.8 


v infaii (km/s) 


-(3.5-0.5) 









" source description from Motte et al. 2003 

b from Beuther et al. 2002, excepted for W43MM1 

c density at 100 AU 

d Log(T) = a Log(r) + j3 



The line profiles observed are also affected by gas motions, 
particularly through the turbulence which is mainly responsible 
for the line emission widths seen in massive protostars. Thus we 
add a turbulent velocity field characterized by a rms velocity vj, 
for each transition of each molecular species (actually OCS lines 
are the most sensitive to the gas components in our model). In 
the case of W43MM1, where asymmetric profiles possibly due 
to infall are detected, we introduce a radial velocity field that 
follows a law in the form u m f = vq/ ^r/r out , and fit the asym- 
metry parameter with different values of vq for each transition 
and molecular species. Concerning the other sources, when the 
line emission profiles are clearly asymmetric, we are not able to 
reproduce those line profiles because: 

- the modelling does not take into account the sub-structure 
(i.e. multiplicity of the massive dense cores, as for 
IRAS05358+3543) which is hence not reproduced correctly 
in the derived source model, 

- outflows are not included in the model while they strongly 
influence the line profiles. 

Some results for CS in the source IRAS 18264-1 152 are shown 
in Fig. 

6.3. Uncertainties 

An overview of uncertainties in the modeling process that we 
used shows that precision on our abundance results are influ- 
enced by four main points: the signal-to-noise ratio, the total 
mass of the source in the model, its temperature and the pop- 
ulations of energy levels of the molecule studied. Furthermore, 
an absolute error comes from the uncertainty on the dust opac- 
ity at millimeter wavelengths, but it disappears in the context of 
our work where the use of the same radiative transfer code and 
the same modeling process keep the comparison between each 
sources relevant. While signal-to-noise ratio is directly known 
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Table 5. Molecular abundances derived for all species within the 4 sources. The temperature of the layer where most of the emission 
comes from is given. The column 4 gives the line whose emission is the strongest at this distance. Abundances for other sources 
(from the literature) are given in the columns 9, 10 and 1 1 for comparison. AFGL2591 is a mid-IR bright HMPO, G31.41+0.31 a 
hot core, and IRAS 16293 a low-mass protostar. 



Source 








W43MM1 


18264 


05358 


18162 


AFGL2591" 


G31.41" 16293"' 


Abundance 


T 


d d 


Line 














(xitr 10 ) 


(K) 


(10 3 AU) 
















*(CS) 


60 


5.9 


3-2 


36 


4.0 


2.3 


- 


100 


70 5100 




75 


4.9 


5-4 


30 


0.26 


1.3 


1.0 








96 


3.3 


7-6 


160 


1.1 


- 


4.6 








60 


5.9 


3-2 


1.8 


0.32 


0.27 


- 








96 


3.3 


7-6 


8.0 


0.50 


- 


0.27 






*(OCS) 


72 


5.2 


8-7 


92 


2.5 


2.4 


1.4 


100 


40 10000 




76 


4.8 


13- 12 


37 


2.4 


0.92 


1.2 








96 


3.3 


19-18 


130 


1.1 
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" Mean abundances from van der Tak et al. 2003 
b Mean abundances from Hatchell et al. 1998 
c Mean abundances from Wakelam et al. 2004a 

d indicative corresponding distance to the center for the source IRAS 18 162-2048 



from observations, the other uncertainties have to be detailed and 
derived to determine how important they are. 

At first, mass uncertainty comes from measurements of the 
flux at millimeter wavelengths. Considering calibration errors, 
we can assume that this value is known with a 20 % accu- 
racy. Error on temperature is directly linked to this first error, 
in addition with the little variation of the bolometric luminosity 
induced. Nevertheless, results of the radiative transfer models 
show that extreme cases make the temperature varying by 3 % 
at the maximum, and tend to compensate first error in the total 
mass. Even if temperature does not change a lot, its influence on 
energy level populations, combined with total mass (hence den- 
sity) uncertainty, can be strong enough to be significant. These 
uncertainties are dependent on the molecule studied and the tran- 
sition observed. As we noted previously, the errors can compen- 
sate themselves, especially when variation of the level energy 
populations of a transition are going the same way. 

Finally, we derived rough values of uncertainties by testing 
temperature and energy level population variations induced by 
flux errors at millimeter range. Results show that our modeling 
process has a relative precision between 25 % and 30 %, depend- 
ing on the considered molecular line emission. This range of 
values does not take into account the signal-to-noise ratio of the 
observations. Concerning SO2 transitions, populations of the en- 
ergy levels are varying widely, making the uncertainty increase 



up to 40 %, even 70 % for the 5 1,5 - 4o i4 where populations 
have an opposite behaviour against temperature variations in our 
models. Furthermore, we note that optically thin emission was 
assumed, whereas some of them are clearly thick (see Table [TOt. 
That is why abundances derived in this case must be treated as 
order of magnitude indicative values. 

6.4. Molecular abundances: results 

The derived abundances for each molecule and each source are 
given in Table[5]for different layers. Each layer is characterized 
by a temperature (and of course by a density, but the population 
of high-energy levels is sensitive to the temperature for dense 
gas tracers). We also give the strongest line emission in that layer 
(but of course lines are excited over several layers). Results are 
plotted in Fig. Inversus the energy of the upper level of the tran- 
sition. Obviously emission lines coming from warmer regions 
(more than 120 K), e.g. involving higher upper energy levels, 
are missing to probe the most inner regions. 

The H2S abundance is roughly the same in the 4 objects 
studied (1.6 - 3.6 10 -10 ) in the outer region (r ~ 60 K), 
showing more variation (0.16 - 4 10~ 10 ) in the warmer inner 
parts (T ~ 100 K): in IRAS 18264-1 152, IRAS05358+3543 
and IRAS18162-2048, the abundance is one order of magnitude 
lower while it remains constant in W43MM1. 
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Fig. 9. Abundances of CS, OCS, SO and S0 2 (relative to H 2 ) 
versus the energy of the upper level of the transition given in 
Table |5]for the 4 sources. 



An obvious trend is seen for SO2, likely because the energy 
ladder of this molecule is best sampled. Its abundance is roughly 
constant (or decreases for W43MM1) in colder layers (T < 80 
K, e.g. > 5 10 3 AU for the source IRAS18162-2048, and E up < 
35 K), then increases by one order of magnitude in warmer lay- 
ers (hence outer to more inner parts), e.g. ~ 3 10 3 AU for the 
same source. Concerning all the other molecules, as fewer tran- 
sitions have been observed, i.e. the energy ladder is not com- 
pletely covered, no obvious trends are seen. In spite of an in- 
sufficient energy coverage, the derived OCS abundance might 
weakly follow the same trend as seen for CS, except in the source 



IRAS18264-1152. We note that the sources IRAS 18264-1 152 
and IRAS05358+3543 exhibit the same abundance variations for 
SO and S0 2 . 

In the three less evolved objects, CS and OCS are the most 
abundant species in the colder outer region while it is SO for 
IRAS 18 162-2048 (for that source OCS is less abundant than SO 
and no CS 3-2 has been observed), SO2 being the least abundant. 

Reproducing the whole CS and C 34 S emission is impos- 
sible with the same isotopic ratio for all the sources. The de- 
rived isotopic ratio [ 32 S/ 34 S] is 9 - 12, 10 - 15, 15 - 20 
and 20 respectively for IRAS05358+3543, IRAS 18264-1 152, 
IRAS 18 162- 2048 and W43MM1 . Th e Solar system abu ndance 
ratio is 22.5 dKahane et al.lll988l) . but IChin et all d 19961) found 
ratios increasing from 14 to 35 with t he distance to the Ga lactic 
center, but this is in contradiction with lNilsson et alj (|2000) who 
derived a ratio of 1 1 towards the very distant source W49N. 

7. Discussion 

7. 1 . W43MM1 : a special case 

W43MM1 is clearly different from the other sources, with the 
highest abundances (relative to H 2 ) of CS, OCS, SO (except in 
the flattening region); the S0 2 abundance in W43MM1 is the 
highest in the warmer inner part. Moreover the plots (see Fig.fTOl 
andfTTTi of the CS and OCS molecular ratios relative to H 2 S un- 
derline this difference, showing high ratios for that source com- 
pared to the others. 

The reason might be that W43MM1 though being a very 
young source has already developed a hot core with t emperatures 
larger than 200 K (e.g. CH3CN lines observed by iMotte et all 
120031) . This hot core region might also be at the origin of the 
unidentified lines described in Sect. 14. II The presence of this hot 
core might also explain the high abundance of H 2 S and OCS 
molecul es which evaporation from grain surface could be im- 
portant dWakelam et al.l2004al) . and consequently the high abun- 
dance of SO. The CS case is different because this molecule is 
already present in the gas before the evaporation. Actually the 
CS chemistry is closely linked to the atomic carbon abundance 
which might be higher in W43MM1 than in the other sources. 

The very high S0 2 abundance close to the center of 
W43MM1 may be due to strong shocks occurring in the inner 
part of the envelope. The fact that W43MM1 is a very young 
source (as indicated by the possibly detected infall, see Sect. 
15.2b suggests that this massive dense core harbors one or sev- 
eral high-mass protostars with powerful outflows, which might 
indu ce shocked regions. S trong shocks can convert SO into SO? 
(see lHatchell et all 1 19981) . leading then to a ratio of unity, as 
observed in the other sources where the process has already 
taken place, transforming the ratio all across the envelope. Hence 
shock chemistry might contribute to the S0 2 formation process. 

Actually, IRAS 05358+3543 harbors a hot core too 
dLeurini et al] |2007|). but very likel y more compact than the 
one of W43MM1 (IMotte et ai]|2003l) . hence more d i luted in the 
beam of our observations. Moreover. iLeurini et al] (|2007) sug- 
gested that the IRAS05358+3543 hot core might not produce a 
very rich chemistry (indeed, these authors underlined differences 
between its molecular spectrum and the one of typical massive 
hot cores). 

7.2. Molecular ratios as a diagnostic tool of evolution ? 

Sulphur is released into the gas phase in the form of H 2 S, which 
is subsequently transformed into SO and later into S0 2 via 
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Fig. 10. Abundance ratio of OCS, CS, SO and S0 2 relative to 
H2S for each source in the layer at T — 60 K. The sources are 
ordered according to the SED classification. 
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Fig. 11. Abundance ratio of OCS, CS, SO and S0 2 relative to 
H2S for each source in the layer at T — 100 K. The sources are 
ordered according to the SED classification. 



neutral-neutral reactio n at 100 K on time scales of ~ 10 3 yr 
dWakelam et al.ll2004al) . The initial destruction of H 2 S by H 3 + 
is even more efficient in the high-mass objects than in low-mass 
proto stars, as water is more abundant here (see lvan der Tak et alj 
120061) . At higher temperature (T> 300 K), H 2 S can be produced 
again. Once in the gas phase, OCS is destroyed later than H 2 S 
(by S + and cosmic rays). CS destruction at low temperature pro- 
duces SO and S0 2 too, while CS is produced again at higher 
temperature from OCS. The time dependence of these process 
might allow us to see an evolution of the abundance along the 
protostar evolutionary track. 

Actually no evolution is obviously seen for the abun- 
dances relative to H 2 , even if they differ significantly from 
source to source. Nevertheless, W43MM1 shows larger abun- 
dances ofCS, OCS and SO compared to IRAS 18264-1 152 and 
IRAS05358+3543 (see Sect. 6.4). 

We have investigated possible correlations of the molecular 
abundance ratios with the order of evolution of our 4 objects for 
two different layers: one around T — 60 - 70 K and a more inner 
one around T — 100 K. To be coherent, we only consider values 
derived for the same layer, i.e same temperature {e.g. transitions 



CS 5 - 4 and OCS 13-12 with E up = 24-37 K, hence mostly 
excited in the 75 K layer). Figures [10] QT| Q~2] and [T3] show that 
some global trend appears, either a decrease or an increase with 
the evolution, but depending on the layer. 

We first compare the molecular abundances relative to H 2 S 
whose abundance is assumed to remain globally constant dur- 
ing the studied evolutionary period and hence comparable from 
source to source. In the colder region (T = 60 K, Fig. ITOb , 
the [OCS/H 2 S] ratio is clearly decreasing from W43MM1 to 
IRAS 18 162-2048 by a factor 10. This trend is also visible for 
CS. The [SO/H 2 S] and [S0 2 /H 2 S] ratios remain roughly con- 
stant; a small increase might actually be present for SO. In 
the inner, hence warmer, part of the envelope (T = 100 K, 
Fig.fTTTi. drawing conclusions is more difficult, because not all 4 
sources were observed in the same transitions. Moreover, again, 
W43MM1 stands out as it shows the highest ratio (except for 
SO), independently of the trend seen for the three other ob- 
jects. Actually, according to the temperature and density distri- 
butions derived by our model (see Fig. |7J, the 100 K layer for 
W43MM1 is closer to the center of the massive core than in the 
other sources, and moreover corresponds to a higher density. But 
the other reason could be that the hot core is dominating. Hence, 
considering the three other sources, an increase of the different 
ratios relative to H 2 S from IRAS 18 1264 to IRAS 18 162-2048 
is seen. Comparing the 60 and 100 K regions, the main differ- 
ence is for OCS and CS: ratios are decreasing at T = 60 K, but 
increasing at higher temperature. 

Relative to OCS (Fig. [T2l . the abundance of CS does not 
show any trend, neither at 76 K nor at 96 K. At 70 - 75 K, 
the [SO/OCS] and [SO9/OCS] ratios increase by at least one or- 
der of magnitude along the derived evolutionary track. At higher 
temperature, hence in the inner regions, the same trend is ob- 
served for the [SO/OCS] ratio, while [S0 2 /OCS] remains ap- 
proximately constant. 

For both probed layers, the [CS/SO] and [S0 2 /SO] ratios de- 
crease (Fig. [TBI, the variation being more important at 60 K. 

Cross-comparison of the plots for outer regions reveals that 
the OCS (and maybe CS) abundance (relative to H 2 S) decreases 
while the SO and S0 2 relative abundance (compared to other 
molecules) increase along the object evolution. This trend is less 
obvious at higher temperature, maybe because i n these regions 
the H 2 S transformation process into SO and S0 2 ( Wakelam et al.l 
l2004alfbl) is becoming less efficient, as it gets close to the comple- 
tion: the new SO and S0 2 molecules have already been produced 
in the inner parts. This hypothesis is strengthened by the larger 
S0 2 abundance measured in all the objects for the larger temper- 
atures, i.e. in warmer environments. Nevertheless, any conclu- 
sions for inner parts are tricky, because of the incomplete energy 
ladder coverage of our observations and so the lack of very high 
energy transitions. 

7.3. Comparison with previous studies 

A comparison between abundances derived for low- and high- 
mass protostars is not straightforward, because of the different 
composition of the ices (from which H 2 S evaporates) and of the 
different physical conditions (density and temperature; e.g. den- 
sity can be two orders of magnitude higher for the same tem- 
perature) as stressed by IWakelam et al.l (l2004a|). Nevertheless, 
Table [5]also gives abundances for AFGL2591 ( van der Tak et al.l 
120031) . a mid-IR bright HMPO (L bo i 2-10 times larger, warmer 
SED with higher F\2 and F^l F tota \), henc e more evolved ob- 
ject, and G3 1 .4 1 + 0.3 1 (iHatchell et al.ll 1 9981) . a massive hot core 
and IRAS 16293 dWakelam et alj|2004al) . We stress that those 
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Fig. 12. Abundance ratio of CS, SO and S0 2 relative to OCS 
for each source in two different layers. The sources are ordered 
according to the sED classification. 



values are mean abundances across the whole envelope. Values 
for AFGL2591 and G31.41 differ a lot for H 2 S and S0 2 (by 
at least a factor 10) but are comparable for CS, OCS, and SO 
to the abundances found in W43MM1. The other sources from 
our sample exhibit abundances one or two order of magnitude 
lower. Abundances in low-mass protostar are 100 - 1000 times 
larger. Clearly the same set of lines has to observed towards these 
sources and the same modelling applied in order to be able to be 
compared with our sample. 

Th e sulphur chemistry was studied by Ivan der Tak et al.l 
d2003l) in the envelopes of more evolved massive protostellar 
object, mid-IR bright HMPOs. The abundances they derived are 
10 - 100 times higher than ours (see Table |5). They concluded 
that shock chemistry is unlikely to contribute and that OCS is a 
major sulphur carrier in the ices. The molecular ratios computed 
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Fig. 13. Abundance ratio of CS and SO2 relative to SO for each 
source in two different layers. The sources are ordered according 
to the sED classification. 



from their work do not show significant variations from source 
to source. 

From our work, the total abundance of sulphur 
(CS+SO+SO2+H2S+OCS) can be estimated to be 
2 10~ 9 - 5 10~ 8 , several orders of magnitude lower than 
the corresponding solar abundance (3.4 10~ 5 ) and the value 
derived in the low-mass hot core IRAS 16293 by Wa kelam et al.l 
(l2004al) (2.8 10~ 6 ). It is likely that we do not probe the very 
inner parts of the protostars where the S-bearing molecules 
are fully evaporated from the grain mantles. Hypothesis to 
explain such depletion of sulphur (including the one previously 
mentioned) will be explored in another paper on the chemistry. 



8. Conclusion 

We have presented here new observations of sulfur-bearing 
species toward 4 massive dense cores. Different gas com- 
ponents, revealing outflows or multiplicity, are detected in 
IRAS05358+3543, IRAS18264-1152 and IRAS 18 162-2048. 
A strong infall of more than 2 kms -1 towards W43MM1 might 
be observed, leading to an impressive kinematic mass infall rate. 
W43MM1 study shows that hot core may appear earlier than ex- 
pected in the evolution. 

An evolutionary classification derived from the calculated 
SEDs is proposed, from W43MM1 to IRAS 18 162-2048, and is 
likely to be observed in the molecular abundances too. More pre- 
cisely, molecular ratios like [OCS/H 2 S], [CS/H 2 S], [SO/OCS], 
[SO2/OCS], [CS/SO] and [SO2/SO], using low energy transi- 
tions, might be good indicators of evolution depending on lay- 
ers probed by the observed molecular transitions: obvious trends 
along the massive core evolution are seen for colder outer re- 
gions at 60 K, but drawing definitive conclusions is more diffi- 
cult for the inner regions. Observations of molecular emission 
from warmer layers, hence implying higher upper energy lev- 
els are mandatory to include. A follow-up with APEX telescope 
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should be made to have access to these higher frequencies. Of 
course, the source multiplicity due the distance is problematic 
and only ALMA will provide a definitive answer. 

Nevertheless, a specific chemical modelling will be done for 
these sources in a forthcoming paper (Wakelam et al., in prepa- 
ration) in order to explain the observed molecular ratios and to 
try to date the studied objects. 
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Fig. 14. Normalized emission of CS, H2S, OCS, SO, SO2 (and isotopic species) lines from W43MM1. Spectra velocity resolutions 
are 0.10-0. 19 km s _1 . The dashed line shows the source LSR velocity. 
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Fig. 15. Normalized emission of CS, H2S, OCS, SO, SO2 (and isotopic species) lines from IRAS18264-1 152. Spectra velocity 
resolutions are 0. 10-0. 19 km s . The dashed line shows the source LSR velocity. 
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Fig. 16. Normalized emission of CS, PbS, OCS, SO, SO2 (and isotopic species) lines from IRAS05358+3543. Spectra velocity 
resolutions are 0. 1 0-0. 1 9 km s~ 1 . The dashed line shows the source LSR velocity. 
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Fig. 17. Normalized emission of CS, H 2 S, OCS, SO, SO2 (and isotopic species) lines from IRAS18162-2048. Spectra velocity 
resolutions are 0. 10-0. 19 km s . The dashed line shows the source LSR velocity. 
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Table 6. Observed line emissions parameters for W43MM1. Widths (Ai^ is the linewidths at 3<x level; Avi/2 is the half 
power linewidth) and line fluxes are derived from gaussian fits (several components) made with CLASS software. For the non 
detection, l<x upper limits are given. 



Species 


Transition 


Av 


T m t> 


<r 


Tt^ — z 

J T mb ou 


A[) 3tr 


Aui/2 


comment 






km/s 


K 


mK 


K.km/s 


km/s 


km/s 




CS 


3-2 


9.3(0.2) 


4.01 


36 


22.90(0.15) 


15.69 


9.21 


a 




5-4 


9.68(0.04) 


6.85 


62 


61.6(0.3) 


14.50 


9.36 


a 




7-6 


4.2(0.3) 


1.2 


255 


5.4(0.7) 


4.5 


4.5 


b 


C J4 S 


3-2 


5.69(0.02) 


2.12 


42 


12.96(0.06) 


11.23 


6.02 






7-6 


6.1(0.3) 


1.05 


269 


4.8(0.2) 


2.7 


5.6 




so 


3 4 -2 3 


6.10(0.01) 


2.66 


23 


17.52(0.03) 


13.24 


6.62 






5 6 -4 5 


4.40(0.06) 


1.14 


65 


5.10(0.07) 


7.22 


4.51 






65 - 5 4 


2.6(0.1) 


0.8 


85 


2.01(0.09) 


3.47 


2.54 






8 8 -7 7 


5.4(0.3) 


0.6 


184 


2.8(0.1) 


5.4 


5.4 


weak det. 


j4 SO 


3 4 -2 3 
8 8 -7 7 


5.6(0.1) 


0.39 


26 
197 


2.42(0.04) 
< 0.009 


10.22 


5.52 


no det. 


so 2 


5 1.5 - 4o,4 


5.05(0.1) 


0.52 


38 


2.33(0.05) 


8.63 


3.91 






lOo.io - 9i >9 
11 1,11 _ lOo.io 


5.1(0.3) 


0.33 


43 
47 


1.41(0.06) 
<0.006 


8.02 


4.0 


no det. 




14o,14 _ 13l,13 






60 


<0.01 






no det. 




18o,18 _ 17l,n 


7.8(0.3) 


0.23 


80 


1.89(0.07) 


7.8 


7.8 






284,24 - 283,25 






60 


<0.01 






no det. 


j4 so 2 


5 1,5 - 4o,4 






32 


<0.006 






no det. 


ocs 


8-7 


6.51(0.02) 


0.90 


11 


6.20(0.02) 


13.10 


6.55 


a 




13 - 12 


4.57(0.03) 


1.58 


39 


7.33(0.05) 


8.72 


4.51 






19-18 


5.37(0.02) 


1.92 


35 


10.46(0.04) 


11.03 


5.01 




OC ,4 S 


8-7 


7.26(0.03) 


0.14 


7 


1.13(0.01) 


13.63 


6.92 




H 2 S 


1 1,0 - lo.i 


9.1(0.2) 


1.06 


72 


12.52(0.6)'' 


12.50 


2.76 


a 




22,0 _ 2[ 1 


4.15(0.04) 


1.26 


45 


5.40(0.05) 


7.62 


4.11 




H 2 J4 S 


1 1,0 - lo.i 






72 


<0.01 






no det. 



" self-absorption at u S0U rce 

b line blend with H 2 CO line from image band. 
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Table 7. Observed line emissions parameters for IRAS18264-1 152. Widths (Avfo is the linewidths at 3<x level; Av\/2 is the 
half power linewidth) and line fluxes are derived from gaussian fits (several components) made with CLASS software. For the 
non detection, l<x upper limits are given. 



Species 


Transition 


Av 


T ra b 


cr 


c 

J T mb <5v 


Al>3o- 


Aui/2 


comment 






km/s 


K 


mK 


K.km/s 


km/s 


km/s 




CS 


3-2) 


4.19(0.01) 


7.2 


54 


31.49(0.07) 


10.28 


4.01 






5-4 


2.97(0.03) 


2.26 


45 


6.81(0.05) 


7.32 


2.81 






7-6 


3.2(0.8)? 


0.6? 


500 


1.9(0.5)? 


? 


? 


b 




3-2 


3.67(0.04) 


1.64 


42 


6.02(0.05) 


8.32 


3.41 






7-6 


5.1(0.5) 


0.34 


111 


1.1(0.1) 


3.7 


4.7 


tentative det. 


so 


3 4 -2 3 


3.27(0.02) 


3.11 


44 


10.38(0.05) 


7.32 


3.21 






5 6 -4 5 


2.98(0.02) 


2.33 


53 


7.08(0.05) 


6.42 


2.81 






6 5 - 5 4 


6.2(0.1) 


1.72 


45 


11.8(0.2) 


11.36 


7.62 


a 




8g - 77 


4.1(0.4) 


0.44 


116 


1.3(0.1) 


2.2 


4.8 


weak det. 


j4 SO 


3 4 -2 3 
8g - 7 7 


3.02(0.2) 


0.19 


39 
132 


0.60(0.04) 
< 0.006 


4.01 


4.01 


d 

no det. 


so 2 


5 1.5 - 4o.4 


3.30(0.3) 


0.24 


33 


0.91(0.05) 


5.91 


5.59 






lOo.io - 9i S 


2.9(0.2) 


0.19 


34 


0.50(0.03) 


3.03 


3.03 


weak line 




11 1,11 ~~ lOo.io 


4.1(0.1) 


0.40 


41 


1.56(0.04) 


7.02 


3.31 






14o,14 _ 13l,13 


3.9(0.3) 


0.34 


59 


0.99(0.07) 


3.95 


3.95 






18o,18 _ 17l,l7 






79 


< 0.004 






no det. 




284,24 - 2 83, 2 5 






63 


<0.01 






no det. 


j4 so 2 


5 1.5 — 4 .4 






64 


<0.01 






no det. 


ocs 


8-7 


3.73(0.1) 


0.18 


5 


0.71(0.01) 


9.19 


3.51 






13 - 12 


3.3(0.1) 


0.41 


46 


1.35(0.05) 


4.13 


3.58 






19-18 


2.7(0.1) 


0.20 


34 


0.49(0.03) 


3.21 


3.21 


bad data 


OC j4 S 


8-7 






4 


<0.001 






no det. 


H 2 S 


1 1,0 - 1(1,1 


3.63(0.04) 


1.94 


59 


7.32(0.06) 


7.22 


3.61 






2 2 ,o _ 2i.i 


2.5(0.3) 


0.14 


28 


0.33(0.03) 


5.51 


2.81 


c 


H 2 J4 S 


1 i,o - lo.i 


2.7(0.1) 


0.35 


60 


0.68(0.06) 


3.72 


3.72 





" self-absorption at u S0U rce 

b line blend with H 2 CO line from image band. 

c unidentified line @ 34.87 km/s (216.716 GHz). 

d + other line @ 29.1 kms" 1 (135.78223 GHz), T mb = 0.14 K, Sv = 1.1 kms 4 , fT mb Sv = 0.17 K-km-s" 1 
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Table 8. Observed line emissions parameters for IRAS05358+3543. Widths (Ai^ is the linewidths at 3<x level; Av\/2 is the 
half power linewidth) and line fluxes are derived from gaussian fits (several components) made with CLASS software. For the 
non detection, l<x upper limits are given. 



Species 


Transition 


Av 


T m b 


a 


7" 

J T mb <5u 


Av 3lT 


AVl/2 


comment 






km/s 


K 


mK 


K.km/s 


km/s 


km/s 




CS 


3-2 


5.54(0.06) 


7.92 


35 


29.75(0.02) 


11.07 


2.77 


a 




5-4 


3.60(0.10) 


5.14 


59 


16.64(0.07) 


7.37 


3.31 


a 




7-6 














no data 


C J4 S 


3-2 
7-6 


4.30(0.16) 


1.59 


24 


4.10(0.01) 


7.12 


2.51 


a 

no data 


so 


3 4 -2 3 


3.97(0.05) 


2.42 


45 


11.46(0.01) 


8.42 


4.71 


a 




5 6 - 4 5 


3.31(0.05) 


2.29 


63 


7.09(0.06) 


6.62 


2.91 


a 




6 5 - 5 4 


2.81(0.08) 


0.99 


60 


2.82(0.06) 


4.54 


2.81 


a 




8g - 7 7 














no data 


,4 SO 


3 4 -2 3 
8 8 -7 7 


4.47(0.35) 


0.26 


50 


1.1 


4.13 


4.4 


a 

no data 


so 2 


5 1.5 - 4o.4 


4.9(0.2) 


0.31 


36 


1.48(0.05) 


5.91 


4.55 






10o,io - 9i, 9 


3.9(0.4) 


0.24 


49 


0.70(0.06) 


1.90 


4.81 






Hi, 11 _ 10o,io 






55 


<0.01 






no det. 




14o,14 - 13l,13 


3.89(0.3) 


0.50 


68 


1.49(0.09) 


4.11 


3.41 






18o,18 _ 17l,i7 














no data 




284.24 - 283, 2 5 






45 


<0.008 






no det. 


,4 so 2 


5 1.5 _ 4 ,4 






25 


<0.005 






no det. 


ocs 


8-7 


4.69(0.05) 


0.27 


10 


1.29(0.01) 


12.9 


5.01 


a 




13- 12 


3.08(0.05) 


0.20 


27 


0.58(0.05) 


3.11 


3.01 


a 




19-18 


5.2(0.2) 


0.46 


58 


2.01(0.06) 


5.61 


5.51 


a 


OC ,4 S 


8-7 






10 


<0.003 






no det. 


H 2 S 


1 1,0 _ lo.i 


4.8(0.15) 


2.56 


53 


9.36(0.15) 


8.46 


3.03 


a 




2 2 ,o _ 2j.i 


2.98(0.15) 


0.37 


44 


0.89(0.04) 


3.91 


3.41 


a 


H 2 j4 S 


1 1,0 _ lo.i 


1.2(0.3) 


0.24 


77 


0.25(0.05) 


0.83 


1.01 


very weak line 6 



" self-absorption at u S0U rce 

* detected at 2cr with two peaks @ -15.57 and -16.21 km/s 
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Table 9. Observed line emissions parameters for IRAS 18 162-2048. Widths(Ai)3 (T is the linewidths at 3<x level; Aui/2 is the 
half power linewidth) and line fluxes are derived from gaussian fits (several components) made with CLASS software. For the 
non detection, l<x upper limits are given. 



Species 


Transition 


Av 


T m b 


<x 


7" 

J T mb ou 


Au 3(r 


An 1/2 


comment 






km/s 


K 


mK 


K.km/s 


km/s 


km/s 




CS 


3-2 














no obs. 




5-4 


2.93(0.01) 


7.41 


38 


22.10(0.03) 


6.78 


2.91 






7-6 


2.46(0.04) 


4.05 


237 


9.3(0.1) 


4.08 


2.14 




C J4 S 


3-2 














no obs. 




7-6 


2.1(0.3) 


0.47 


116 


0.87(0.06) 


1.17 


2.63 




so 


3 4 -2 3 


2.24(0.01) 


5.27 


40 


12.00(0.04) 


6.46 


2.10 






5 6 - 4 5 


2.20(0.01) 


5.22 


45 


11.31(0.03) 


5.42 


4.84 






6 5 - 5 4 


2.24(0.03) 


2.80 


79 


6.25(0.08) 


4.52 


1.78 






8g - 7 7 


2.2(0.2) 


1.04 


230 


1.7(0.1) 


1.1 


1.7 




"SO 


3 4 -2 3 
88 - 7 7 






38 
123 


<0.007 
< 0.006 






no det. 
no det. 


so 2 


5 1.5 - 4o,4 


3.3(0.1) 


0.72 


34 


1.86(0.04) 


6.46 


1.68 






10 ,,o - 9i, 9 


1.69(0.07) 


0.43 


33 


0.75(0.02) 


2.39 


1.75 


a 




II1.11 - 10o,io) 


3.86(0.08) 


1.14 


35 


4.20(0.07) 


7.88 


3.49 






14o,14 - 13 1,13 


3.63(0.09) 


1.02 


59 


3.71(0.07) 


5.81 


3.68 






18o,18 ~~ 1 '1,17 


1.2(0.5) 


0.13 


75 


0.16(0.05) 




3.7? 


tentative det. 




28 4 . 24 - 283,25 






42 


<0.007 






no det. 


J4 so 2 


5 1,5 _ 4 ,4 






38 


<0.007 






no det. 


ocs 


8-7 


3.5(0.1) 


0.12 


8 


0.40(0.01) 


7.17 


3.10 






13 - 12 


4.7(0.5) 


0.23 


51 


0.75(0.06) 


2.85 


6.26 






19-18 


5.12(0.08) 


0.55 


33 


2.53(0.04) 


7.55 


4.55 




OC J4 S 


8-7 






9 


<0.003 






no det. 


H 2 S 


1 1,0 - lo.i 


3.46(0.02) 


2.92 


50 


10.00(0.05) 


8.0 


3.36 


a 




2 2 ,o _ 2i,i 


3.2(0.3) 


0.24 


24 


0.78(0.07) 


4.26 


2.97 


b 


H 2 J4 S 


1 1,0 _ lo.i 


2.3(0.1) 


0.60 


65 


1.15(0.06) 


2.58 


1.55 





" self-absorption at u S0U rce 

b unidentified lines @ 216.701026 and 216.71.6 GHz (25 and 4 km/s) 



Table 10. Opacities derived from the modelling of the molecular emission of the sources. "ND" stands for no data. 



Lines 


W43-MM1 


18264 


05358 


18162 


CS 3-2 


52.0 


1.6 


0.45 


ND 


CS5-4 


65.0 


0.31 


0.81 


0.67 


CS 7-6 


300 


0.94 


ND 


2.8 


C J4 S 3-2 


2.8 


0.18 


0.09 


ND 


C 34 S 7-6 


16.0 


0.51 


ND 


0.41 


OCS 8-7 


6.0 


0.07 


0.04 


0.06 


OCS 13-12 


2.6 


0.1 


0.06 


0.1 


OCS 19-18 


4.9 


0.12 


0.07 


0.2 


OC ,4 S 8-7 


0.4 


0.03 


0.02 


0.05 


H2S 1 1.0 - lo.i 


4.2 


0.53 


0.21 


3.9 


H2S 2 2 .o - 2|.i 


3.3 


0.13 


0.08 


0.15 


H 2 j4 S 1 1.0 - lo.i 


0.36 


0.11 


0.04 


0.48 


SO 3 4 - 2 3 


6.3 


0.25 


0.32 


0.36 


SO 5 6 - 4 5 


0.72 


0.21 


0.15 


0.39 


SO 6, - 5 4 


0.52 


0.21 


0.10 


0.28 


SO 8 8 - 7 7 


6.8 


0.33 


ND 


0.59 


J4 SO 3 4 - 2 3 


0.45 


0.05 


0.11 


0.89 


34 SO 8 8 - 7 7 


1.6 


0.26 


ND 


2.9 


SO, 5,,5 - 4 ,4 


0.42 


0.06 


0.03 


0.08 


SO, 10 ,io - 9i, 9 


0.33 


0.06 


0.03 


0.04 


S0 2 ll U i - lOo.io 


ND 


0.11 


ND 


0.13 


SO, 14 ,,4- 13 u3 


ND 


0.21 


0.23 


0.42 


S0 2 18o. 18 - 17i.i 7 


6.0 


0.24 


ND 


0.22 


J4 S0 2 5 li5 - 4 ,4 


0.18 


0.61 


0.02 


0.03 



